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ABSTRACT 

Identifying the processes that determine strength, duration and variability of protostellar mass growth is a fun- 
damental ingredient of any theory of star formation. I discuss protostellar mass accretion rates M from numerical 
models which follow molecular cloud evolution from turbulent fragmentation towards the formation of stellar 
clusters. In a dense cluster environment, M is strongly time varying and influenced by the mutual interaction of 
protostellar cores and their competition for accretion from the common cluster gas reservoir. Even for protostars 
with similar final mass, the accretion histories may differ dramatically. High-mass stars build up in the central 
parts of dense, cluster-forming cloud regions. They begin to form in the early phases of cluster evolution, and 
continue to grow at a high rate until the available gas is exhausted or expelled by feedback. Lower-mass stars tend 
to form at later phases, and M declines rapidly after a short initial phase of strong growth. I present a simple fit 
formula for the time evolution of the average M for protostars of different masses in a dense cluster environment. 

Subject headings: hydrodynamics - ISM: kinematics and dynamics - stars: formation - stars: pre-main sequence 



1. INTRODUCTION 

Star formation is a complex process. In the standard theory 
of star formation, stars build up from the "inside-out" collapse 
of singular isothermal spheres, which are generally assumed 
to result from the quasistatic contraction of magnetically sup- 
ported cloud cores due to ambipolar diffusion (Shu 1977, Shu, 
Adams, & Lizano 1987). This picture is able to reproduce many 
observed features of protostars, however, has always been chal- 
lenged by more dynamical points of view (see e.g. Whitworth et 
al. 1996 for a critical discussion). An alternative approach takes 
interstellar turbulence into account, where supersonic random 
motions in molecular clouds lead to transient shock-generated 
density fluctuations. Some may exceed the threshold of grav- 
itational stability and collapse to form stars (e.g. Elmegreen 
1993, Padoan 1995, Padoan & Nordlund 1999). The location, 
timescale and efficiency of star formation hereby intimately de- 
pend on the stochastic properties of the underlying turbulent ve- 
locity field (Klessen, Heitsch, & Mac Low 2000, Heitsch, Mac 
Low, & Klessen 2001, hereafter KHM and HMK, respectively, 
also Klessen 2000). The majority of stars form in a clustered 
environment (see the review by Clarke, Bonnell, & Hillenbrand 

2000, and Elmegreen et al. 2000). For massive clusters, the mu- 
tual interaction of collapsing protostellar cores in the dense cen- 
tral region and their competition for accretion from a common 
gas reservoir are important processes which determine stellar 
mass growth (Bonnell et al. 1997, Klessen & Burkert 2000, 

2001, hereafter, KB1 and KB2, respectively). It is this star 
forming environment which is considered in the present paper. 

A gravitationally unstable, collapsing gas clump follows an 
observationally well determined sequence. Prior to the forma- 
tion of a hydrostatic nucleus, an observed pre-stellar condensa- 
tion exhibits a density structure which has a flat inner part, then 
decreases outward roughly as p oc r~ 2 , and is truncated at some 
finite radius (e.g. Bacmann et al. 2000). Once the central YSO 
builds up, the class phase is reached and the density follows 
p oc r" 2 down to the observational resolution limit. As larger 
and larger portions of the infalling envelope get accreted the 
protostar is identified as class I object, and when accretion fades 
away it enters the T Tauri phase (e.g. Andre, Ward-Thompson, 



& Barsony 2000). IR and sub-mm observations suggest that the 
accretion rate M varies strongly and declines with time. Accre- 
tion is largest in the class phase and drops significantly in 
the subsequent evolution (e.g. Andre & Montmerle 1994, Bon- 
temps et al. 1996, Hendriksen, Andre, & Bontemps 1997, here- 
after HAB). The estimated lifetimes are a few 10 4 years for the 
class and a few 10 s years for the class I phase. 

These observational findings favor the dynamical picture 
of star formation (e.g. Larson 1969, Penston 1969, Hunter 
1977, HAB, Basu 1997), but could in principle be recon- 
ciled with an "inside-out" collapse model when adopting ini- 
tial conditions differnt from the classical quasi-static singular 
isothermal sphere. All analytical studies and most numeri- 
cal work of protostellar core collapse (e.g. Foster & Chevalier 
1993, Tomisaka 1996, Ogino, Tomisaka, & Nakamura 1999, 
Wuchterl & Tscharnuter 200 1 , herafter WT) concentrate on iso- 
lated objects. It is thus the aim of the present study to illustrate 
the effect of a dense cluster environment on protostellar mass 
accretion rates. 

2. NUMERICAL MODELS OF CLUSTERED STAR FORMATION 

In the approach adopted here, molecular clouds form clusters 
of stars in regions where interstellar turbulence becomes too 
weak to supply sufficient support against gravitational collapse, 
either because the region is compressed by a large-scale shock 
(KHM, HMK), or because interstellar turbulence is not replen- 
ished and decays on short timescales (Mac Low et al. 1998, 
Stone, Ostriker, & Gammie 1998, Mac Low 1999). The lat- 
ter scenario is investigated by Klessen, Burkert, & Bate (1998), 
KB1, and KB2, using smoothed particle hydrodynamics (SPH) 
simulations of the evolution and fragmentation of molecular 
cloud regions where decaying turbulence is assumed to have 
left behind random Gaussian density fluctuations. The com- 
plete lack of turbulent support leads to rapid star formation in 
a strongly clustered mode. The interplay between self-gravity 
and isothermal gas pressure considered in the models repro- 
duces the basic features of young star clusters remarkably well, 
e.g. the mass distribution of protostars is well fit by a log- 
normal with width and peak comparable to observational es- 
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timates of the IMF (KB1, KB2). The current analysis is based 
on the high-resolution simulation X discussed in KB 1 . 

The considered molecular cloud region has a mean density of 
n(H 2 ) = 10 5 cm -3 and a temperature of 10K. Its initial Gaussian 
density field follows a power spectrum P(k) oc k~ 2 in a volume 
of (0.32pc) 3 containing roughly 200 M Q . As the system con- 
tracts gravitationally, a cluster of 56 protostellar cores builds 
up. As no feedback is included, the common cluster gas reser- 
voir from which these cores feed is exhausted after roughly 
4 x 10 5 years, corresponding to <~ 3 global free-fall timescales. 

Once the nucleus of a collapsing gas clump in the model ex- 
ceeds a density n(H2) ~ 10 9 cm" 3 , it is identified as protostel- 
lar core and replaced by a 'sink' particle, which has the ability 
to accrete matter while keeping track of mass, linear and an- 
gular momentum (Bate, Bonnell, & Price 1994). Matter that 
accretes through the boundary of the 'sink' volume would con- 
tinue to fall and reach the central hydrostatic YSO within less 
than ~ lOOOyears (WT). The diameter of the accretion volume 
is 320 AU, roughly equivalent to the Jeans scale at the threshold 
density. For the stellar mass range considered here feedback 
effects are not strong enough to halt or delay accretion onto the 
stellar photosphere (Wuchterl & Klessen 2001, hereafter WK). 
Thus, the core accretion rates derived here are good estimates 
for the actual stellar accretion rates. Deviations are expected 
only for very high-mass stars, which will at some stage begin 
to emit UV photons, or for protostellar cores forming binary 
stars, where the infalling mass must be distributed between two 
stars, or if very high-angular momentum material is accreted, 
where a certain mass fraction may end up in a circumbinary 
disk and not accrete onto a star at all. 

3. TIME- VARYING PROTOSTELLAR ACCRETION RATES 

The following conclusions about the mass accretion rate in 
dense clusters can be derived: 

(1) Protostellar accretion rates in a dense cluster environ- 
ment are strongly time variable. This is illustrated in Figure 1 
for 49 randomly selected cores. 

(2) The typical density profiles of gas clumps that give birth 
to protostars exhibit a flat inner core, followed by a density fall- 
off p oc r~ 2 , and are truncated at some finite radius (in the dense 
centers of clusters often due to tidal interaction with neighbor- 
ing cores). Figure 12 in KB1 presents examples. As result, a 
short-lived initial phase of strong accretion occurs when the flat 
inner part of the pre-stellar clump collapses. This corresponds 
to the class phase of protostellar evolution. If these cores 
were to remain isolated and unperturbed, the mass growth rate 
would gradually decline in time as the outer envelope accretes 
onto the center. This is the class I phase. Once the truncation ra- 
dius is reached, accretion fades and the object enters the class II 
phase. This behavior is expected from analytical models (e.g. 
HAB) and agrees with other numerical studies (e.g. Foster & 
Chevalier 1993, WT). However, collapse does not start from 
rest for the density fluctuations considered here, and the accre- 
tion rates exceed the theoretically predicted values even for the 
most isolated objects in the simulation. 

(3) The mass accretion rates of cores in a dense cluster de- 
viate strongly from the rates of isolated cores. This is a direct 
result of the mutual dynamical interaction and competition be- 
tween protostellar cores. While gas clumps collapse to build 
up protostars, they may merge as they follow the flow pattern 
towards the cluster potential minimum. The timescales for both 
processes are comparable. The density and velocity structure 



of merged gas clumps generally differs significantly from their 
progenitor clumps, and the predictions for isolated cores are no 
longer valid. More importantly, these new larger clumps con- 
tain multiple protostellar cores, which subsequently compete 
with each other for the accretion from a common gas reservoir. 
The most massive core in a clump is hereby able to accrete more 
matter than its competitors (Bonnell et al. 1997, KB1). Its ac- 
cretion rate is enhanced through the clump merger, whereas the 
accretion rate of low-mass cores typically decreases. Tempo- 
rary accretion peaks in the wake of clump mergers are visible 
in abundance in Figure 1. Furthermore, the small aggregates 
of cores that build up are dynamically unstable and low-mass 
cores may be ejected. As they leave the high-density environ- 
ment, accretion terminates and their final mass is reached. 

(4) The most massive protostars begin to form first and con- 
tinue to accrete at high rate throughout the entire cluster evo- 
lution. As the most massive gas clumps tend to have the largest 
density contrast, they are the first to collapse and constitute the 
center of the nascent cluster. These protostellar cores are fed at 
high rate and gain mass very quickly. As their parental clumps 
merge with others, more gas is fed into their 'sphere of influ- 
ence'. They are able to maintain or even increase the accretion 
rate when competing with lower-mass cores (e.g. core 1 and 8 in 
Figure 1). Low-mass stars, on average, tend to form somewhat 
later in the dynamical evolution of the system (as indicated by 
the absolute formation times in Figure 1; also Figure 8 in KB1), 
and typically have only short periods of high accretion. 

(5) Individual cores in a cluster environment form and evolve 
through a sequence of highly probabilistic events, therefore, 
their accretion histories differ even if they accumulate the same 
final mass. Accretion rates for protostars of certain mass can 
only be determined in a statistical sense. For the investigated 
cluster, I define four mass bins each containing 14 cores, and 
calculate the average accretion rate (M) and its mean absolute 
deviation. The result is shown in Figure 2, with the considered 
mass range indicated at the top of each plot. An exponential 
decline with a cut-off at f en( j offers a reasonable approximation 
to (M), 

log 10 (M) [Moyr 1 ] =A J^EL for < f-f form < f end (1) 

with normalization A and decline time r obtained from fitting 
(M) over 80% of the average accretion duration. There is an im- 
plicit upper mass limit M x = lO'V/lnlO M Q when f en d — ► oo. 
The integration of (M) yields (M en d) < M^. The values A, 
t, (M en d), and Moo are listed in Table 1. Note, as individual 
accretion rates deviate considerably from the cluster average, 
when constructing cluster accretion rates from the fit formula 
a time-dependent random component with roughly 50% devia- 
tion should be superimposed onto the mean. 



Table 1 



A (Moyr" 1 ) 


r(10 5 yr) 


(M end ) (Mq) 


Moo (Mq) 


bin 


-4.13 


2.24 


4.17 


7.21 


(a) 


-4.16 


1.03 


2.04 


3.09 


(b) 


-4.35 


0.79 


1.08 


1.54 


(c) 


-4.66 


0.49 


0.39 


0.47 


(d) 



Fit parameters to (M) applying Equation 1 to the initial 80% 
of the mean accretion duration in mass bins (a) to (d). 
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Figure 1: Examples of time varying mass accretion rates for protostellar cores forming in a dense cluster environment. The left panel shows accretion rate M 
versus time after formation t — tf or m for 49 randomly selected protostellar cores in the high-resolution model I of KB1. The time of core formation ti otm is identified 
by the occurance of the first hydrostatic stellar object. Accretion rates at times t — tf otxa < indicate the mass flow into the control volume prior to that (the gas 
clump is still in the pre-stellar phase). To link individual accretion histories to the overall cluster evolution, tf or m is indicated in the upper right corner of each 
plot and measures the elapsed time since the start of the simulation. The free-fall timescale of the considered molecular region is « 10 5 years. High-mass stars 
tend to form early in the dynamical evolution and are able to maintain high accretion rates throughout the entire simulation. On the contrary, low-mass stars 
tend to form later in the cluster evolution and M declines strongly after the short initial peak accretion phase. Altogether, the accretion histories of cores (even of 
those with similar masses) differ dramatically from each other due to the stochastic influence of the cluster environment, as clumps merge and protostellar cores 
compete for accretion from a common gaseous environment. The right panel plots for the same cores M as function of the accreted mass M with respect to the 
final mass M en d, which is indicated in the center of each plot. Note that the mass range spans two orders of magnitude. 
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Figure 1: Averaged mass accretion rate (M) (thick line) as function of time relative to core formation £ — i form for 
four different mass bins (as denoted in the top section of each plot) overlaid on the contributing individual accretion 
histories. The mean absolute deviation from (M) is indicated by thin lines. The exponential approximation 
(Equation 1) to (M) and the adopted fit interval (i.e. the first 80% of the mean accretion duration) are indicated 
by the dashed line. The four mass intervals contain equal numbers of protostellar cores. 



(6) As M and M en d of evolving protostars in dense clusters 
are influenced by mutual stochastic interactions, the bolomet- 
ric luminosity Lboi and temperature 7boi are not functions of 
mass and age alone (Myers et al. 1998, WT), but also depend 
on the statistical properties of the parental cluster (WK). For 
protostars in the main accretion phase, mass and age determi- 
nation from comparison with theoretical tracks, either in the 
log 7b i — logLbol or the Hertzsprung-Russell (H-R) diagram, is 
only possible in an approximated way - as the average over 
many different theoretical accretion histories for different clus- 



ter environments. This may still affect protostars when they 
enter the pre-main sequence contraction phase and accretion 
fades off. A quantitative analysis of the expected errors for dif- 
ferent cluster environments will be the subject of a subsequent 
investigation. 
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